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Abstract 

The late-time evolution of the charged massive Dirac fields in the background of a Reissner-Norstrom 
(RN) black hole is studied. It is found that the intermediate late-time behavior is dominated by an inverse 
power-law decaying tail without any oscillation in which the dumping exponent depends not only on the 
multiple number of the wave mode but also on the field parameters. It is also found that, at very late times, 
the oscillatory tail has the decay rate of t~^^^ and the oscillation of the tail has the period 27r//x which is 
modulated by two types of long-term phase shifts. 
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I. INTRODUCTION 

The evolution of field perturbation around a black hole consists roughly of three stages The 
first one is an initial wave burst coming directly from the source of perturbation. The second one 
involves the damped oscillations called the quasinormal modes. And the last one is a power-law tail 
behavior of the waves at very late time. 

The late-time evolution of various field perturbations outside a black hole has important impli- 
cations for two major aspects of black-hole physics: the no-hair theorem and the mass-inflation 
scenario 3] There: ore, the decay rate of the various fields ^as been extensively studied j^-js^ 
since Wheeler 0, 24 1 introduced the no-hair theorem. Price |4] studied the massless external per- 
turbations and found that the late-time behavior for a fixed r is dominated by the factor t-(2'+3)_ 
Barack, Ori and Hod considered the late-time tail for the gravitational, electromagnetic, neu- 

trino and scalar fields in the Kerr spacetime. Starobinskii and Novikov P] analyzed the evolution 
of a massive scalar field in the RN background and found that there are poles in the complex plane 

n 

which are closer to the real axis than in the massless case. Hod and Piran 9] pointed out that, if the 
field mass /i is small, the oscillatory inverse power-law behavior $ ~ i~('+3/2) sin(^t) dominates as 
the intermediate late-time tails in the RN background. We 10| studied the late-time tail behavior 
of massive Dirac fields in the Schwarzschild black-hole geometry and found that this asymptotic be- 
havior is dominated by a decaying tail without any oscillation. Koyama and Tomimatsu ll| found 
that the very late-time tail of the massive scalar field in the Schwarzschild and RN background is 
approximately given by t~^/Sin(/it). 

Although much attention has been paid to the investigations of the late-time behaviors of the 
neutral scalar, gravitational and electromagnetic fields in the static and stationary spacetimes, at 
the moment the study of the late-time tail evolution of the charged massive fields is still an open 
question. The main purpose of this paper is to study the late-time tail evolution of the charged 
massive Dirac fields in the RN black-hole background. 



II. LATE-TIME TAIL OF THE CHARGED MASSIVE DIRAC FIELDS 



In the RN spacetime the Dirac equations coupled to a electromagnetic fields j25[ can be separated 
by using the Newman-Penrose formalism 2^ . After the tedious calculation, we find that the angular 
equation is the same as in the Schwarzschild black hole and the radial equation can be expressed 
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where A = r"^ — 2Mr+Q'^ (M and Q represent the mass and charge of the black hole), (ir* = {r^/A)dr 
and \E'± = rA^3(A^ + ^2^2^-ig±iarctan(^)/2^^^ (-^±| an usual radial wave function 0). 

Let us assume that both the observer and the initial data are situated far away from the black 
hole. Then, we can expand Eq. 1)2.11) as a power series in M/r and Q/r and obtain (neglecting terms 
of order O ((cu/r)^) and higher) 
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where w = a = {Mfi^ + eQuj)/uj - 2Mw, = \ + X"^ + AMjj^ - Q^{e^ + /i^) - 2iseQ + 

(X/ H + 8MeQ + 2isM)uj and ^± = (A/r^)^/^\E'-|-. We obtain the two basic solutions required to build 
Green's function 



^1 



Ae-^''{2wr)^+^M{- + b-a,l + 2b, 2wt) 
2 = 5e-^'^(2ti7r)^+^f/(l + 6 - a, 1 + 26, 2vDr), 



(2.3) 



where A and B are normalization constants, M(a, 6, z) and [/(a, 6, 2;) represent the two standard 
solutions to the confluent hypergeometric equation 1271 . U{d,b,z) is a many- valued function, i.e., 
there is a cut in '^2- Hod, Piran and Leaver P, bsl ] found that the asymptotic massive tail is 
associated with the existence of a branch cut (in "$2) placed along the interval —fi < 00 < fi and the 
branch cut contribution to the Green's function is 



G^(r*,r:;t) 



1 
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F(w)e'''^^duj 
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with 



^i«,roe^^)^2(r-.,t^e^ 
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where W{w) = W^^i, ^2) = ^i^2,x - ^2'^i,x is the Wronskian. We obtain, with the help of Eq. 



; 13. 1.22) of Ref. 



I, that W{we™) = —W{'cj) = AB r(i/2+b-a) 4&07. For simphcity we assume that 
the initial data has a considerable support only for r values which are smaller than the observer's 
location. This, of course, does not change the late-time behavior. Noting that when t is large, 
the term e"*"^* oscillates rapidly. This leads to a mutual cancellation between the positive and 
the negative parts of the integrand, so that the effective contribution to the integral arises from 



0{fi - i) or equivalently w = 0{^/Ji/i) Q. Using Eqs. (13.1.32), (13.1.33) and (13.1.34) of 



Ref. (22!, we find that F{zu) is given by 



F{w) = — — - — ^ — -[M{- + 6 + a, 1 + 26, 2ti7<)M(- - 6 + a, 1 - 26, 2^^) - M{- + 6 - a, 1 + 26, 2wr'^) 

xM(i - 6 - a. 1 - 26, 2^r.)] _ ^^'^^^^^^ + ^ ' [M(l + 6 - M + 26, 2^r:) 

4 ' ' r(26)r(i-6-a) 4w6 ^4 ' *' 

xM(i + 6 - a, 1 + 26, 2wr^) + e(i+2^)"A-/(i + 6 + a, 1 + 26, 2vjrl)M{^ + 6 + a, 1 + 26, 2wr^)]. (2.5) 

We first focus our attention on the intermediate asymptotic behavior of the charged massive 
Dirac fields. That is the tail in the range M ^ r ^ t ^ M/(Myu)^. In this time scale, we find 
that the frequency range = 0{^y fi/t), which gives the dominant contribution to the integral, 
implies a ^ 1. Equation ()2.2|) shows that a originates from the 1/r term which describes the effect 
of backscattering off the spacetime curvature. That is to say, the backscattering off the spacetime 
curvature from the asymptotically far regions is negligible for the case a ^ 1. Then, using the fact 
that M(a, 6, 2;) ~ 1 as z — > 0, we have 

sin{Txh) 2i+2''62 r(6)2 
Substituting Eq. (jSIEI) into the Eq. (j231), we find 



= .J_.. ^iry:)y^^\ (2.6) 



G^(n,r:;t) = 

Unfortunately, the integral ()2.7p can not be evaluated analytically since the parameter h depends on 
uj. However, we can work out the integral numerically and the corresponding results are presented in 
Figs. EEl Figure n describes In |G'^(r^,, r'^] t) \ versus t for different eQ, which shows that the dumping 
exponent depends on seQ, i.e., the product of the spin weight of the Dirac fields and the charges of 
the black hole and Dirac fields, and seQ < speeds up the decay of the perturbation but seQ > 
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FIG. 1: Graphsofln|G^(r,,r:^;t)| versus t for different seQ. These figures show that the dumping exponent 
depends on the product of the spin weight and the charges of the black hole and Dirac fields, and seQ < 
speeds up the perturbation decay but seQ > slows it down. 

slows it down. Figure |21 illustrates In \G'-^ {r^,r'^]t)\ versus t for different A with s = ±1/2, fi = 0.01 
and eQ = 0.01, which indicates that the dumping exponent depends on the multiple number of the 
wave mode, and the larger the magnitude of the multiple number, the more quickly the perturbation 
decays. Figure El gives In |G*"(r*, r^; t)| versus t for different mass fi of the Dirac fields, which shows 
that the dumping exponent depends on the mass of the Dirac fields, and the smaller the mass /i, the 
faster the perturbation decays. 

In the above discussion we have used the approximation of a <C 1, which only holds when 
fit < l/ifiM)^. Therefore, the power-law tail found earlier is not the final one, and a change 
to a different pattern of decay is expected when a is not negligibly small. Here we examine the 
asymptotic tail behavior at very late times such that fit ^ l/(/iM)^. This asymptotic tail behavior 
is caused by a resonance backscattering due to spacetime curvature lll|. In this case, we have 
a ^ (MyU^ + eQfi)/zu ^ 1, namely, the backscattering off the spacetime curvature in asymptotically 
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FIG. 2: Graphs of In |G'^(r*, r^; t)\ versus t for different A with s = ±1/2, fi = 0.01 and eQ = 0.01, showing 
that the dumping exponent depends on the multiple number of the wave mode, and the larger the magnitude 
of multiple number, the more quickly the perturbation decays. 

far regions is important. Using Eq. (13.5.13) of Ref. and Eq. (j2.5p . we obtain 

r(l + 26)r(l-26)<r. 



2b 



[J2b{\/or{l)J-2b{.^/^*) - /26(A/m^)/-2fe(yan)] 
r(l + 2b)^r(-2&Kr, V{\ + h-a) 

2hV{2h) 2hT{\-h-af 

[J2biy^,)J2biV^) + ^2fe(v^)^26(v^)], (2-8) 

where a = ^{Mji^ + eQjj), and I±{2b) is the modified Bessel function. 

We can study late-time behavior for the first term of Eq. (j2.8p using numerical method. The 
result is presented by Fig. 0] which shows that asymptotically late-time tail arising from the first 
term is still ~ although the factor F(l + 26)F(1 — 2b) /2b is not a constant. 

Now let us to find the behavior of the second term. Because of F(l + 26)F(1 — 2b) /2b = — F(l + 
26)2F(-26)/(26F(26)) and the factor F(l + 26)F(1 - 26) /26 almost dose not affect the asymptotical 
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FIG. 3: Graphs of In {r^,r'^;t)\ versus t for different mass fi. The result obtain from these figures is that 
the dumping exponent depends on the mass of the Dirac fields, and the smaller the mass n, the faster the 
perturbation decays. 

behavior of the first term, we can define 

r(i + 2fe)^r(-2&Kr, ^ 
" 2bT{2b) [Mv^^OMv^) 

which approximates to a constant. Then, in the limit a ^ 1, the contribution of the second term to 
the Green's function can be expressed as 

G^(r„ <; t) ~ — /"^ e*(2™— *)e^'^(ia;, (2.9) 



27r 



where the phase (j) is determined by e*''' = ^^(^_lyt^2i^a , and it remains in the range < < 27r, 
even if a becomes very large. We use saddle-point integration to evaluate Eq. ()2.9p and find 

G^{n,r:;t) ~ ^ (2vr)t (M/i + 2eQ)^(/it)-t 

sm{/xt - [27r(M/i + 2eg)]3(/xt)3 - ^(cuo) - j}, (2.10) 
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FIG. 4: Graphs of \n\G'^ {r^,r'^;t)\ versus t with M = 100, /x = 0.01, s = -1/2, A = 1, Q = O.IM and 
e = 0.002 for the first term (solid line). The dashed line is ~ log \ l/t\. 



which is the asymptotic behavior of the Green's function at very late times. Eq. ()2.10|) shows that 
the decay rate of the asymptotic tail is t~^^^ and the oscillation of the tail has the period 27r//i 
which is modulated by two types of long-term phase shifts, a monotonously increasing phase shift 

2 1 

[27r (M/i + 2e(5)] 3 (/it)3 and a period phase shift 4>{uJo). 

To confirm the analytical prediction, we present the numerical result of the second term in Fig. 
Eland find that the decay rate of the asymptotic tail is still t~^^^. 



III. SUMMARY 



The intermediate late-time tail and the asymptotic tail behavior of the charged massive Dirac 
fields in the background of the RN black hole are studied. The results of the intermediate late-time 
tail are presented by figures because we can not obtain analytically Green's function G'"{r^,,rl]t) 
because the parameter b in the integrand of the Green's function depends on the integral variable 
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FIG. 5: Graphs of ln\G^ {r^,ri;t)\ versus t with M = 100, /i = 0.01, s = -1/2, A = 1, Q = O.IM and 
e = 0.002 for the second term (solid hne). The dashed hne is ~ log 

cu. We learn from the figures that the intermediate late-time behavior is dominated by an inverse 
power-law decaying tail without any oscillation, which is different from the oscillatory decaying tails 
of the scalar fields. It is interesting to note that the dumping exponent depends not only on the 
multiple number of the wave mode but also on the mass of the Dirac fields and the product seQ. 
We also find that the decay rate of the asymptotically late-time tail is t~^/^ and the oscillation of 
the tail has the period 27r/// which is modulated by two types of long-term phase shifts. 
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